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£\J , Abstract 

^ , The statistical bootstrap model is critically revised in order to include a 

O ■ medium-dependent resonance width in it. We show that a thermodynamic 

' model with a vanishing width below the Hagedorn temperature Th and a 

Hagedorn spectrum-like width above Th may not only eliminate the diver- 
gence of the thermodynamic functions above Th, but also gives a satisfactory 
description the lattice quantum chromodynamics (QCD) data on the energy 
density above the chiral/deconfinement transition as the main result of this 
contribution. This model allows to explain the absence of heavy resonance 
\ contributions in the fit of the experimentally measured particle ratios at SPS 

f^**) . and RHIC energies. 
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1. Introduction 



O 

^ ' Ultrarelativistic heavy-ion collision experiments at SPS and RHIC are performed 

with the aim to create conditions of temperature and density under which hadronic 
matter undergoes a phase transition to the hypothetical quark-gluon plasma (QGP) 
and to investigate the properties of this state of matter once it is created. The 
strongest theoretical support for the existence of the QGP comes from lattice QCD 
simulations at finite temperature T which show a step-like enhancement of the 
effective number of degrees of freedom (energy density in units of T 4 ) at a critical 
temperature T c , see Fig. ^ This behavior is conventionally interpreted as the 
transition from a few mesonic degrees of freedom (mainly 7r and p mesons) to those 
of quasifree quarks and gluons. A microscopic description of this transition is, 
however, still missing. On the one hand, it has been shown that a resonance gas 
model can perfectly explain the steep rise in the number of degrees of freedom at 
T«T C . On the other hand, lattice QCD has also revealed that hadronic correlations 
persist for T > T c 0. The question arises whether it is more appropriate to describe 
hot QCD matter in terms of hadronic correlations rather than in terms of quarks and 
gluons. In the present contribution, we introduce a generalization of the Hagedorn 
resonance gas (statistical bootstrap) model as such a description. 
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Figure 1: Energy density of lattice QCD (circles) and Hagedorn gas (long dashed 
line) as functions of the reduced temperature T/T c . The contributions of pion 
gas (dotted line) and gas of pions and p mesons (dashed line) are shown for a 
comparison. 

The statistical bootstrap model (SBM) 01 El EH is based on the hypothesis 
that hadrons are made of hadrons, with constituent and compound hadrons being 
treated on the same footing. This implies an exponentially growing form of the 
hadronic mass spectrum ph{iti) w Ch^ 11 exp [tti/Th] for m — * oo. The parameter 
Th, Hagedorn temperature, was interpreted as a limiting temperature reached at 
infinite energy density. 

The extensive investigation of the SBM has led to a formulation of both the 
important physical ideas and the mathematical methods for modern statistical me- 
chanics of strongly interacting matter. Thus it has been clarified that the original 
SBM requires two crucial modifications: 

• hadrons should be considered as composite objects (the simplest way are 
hadrons as MIT bags [3| of quarks and gluons) and, their proper volume 
v has to be taken into account [HI E] in the partition function; 
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• only SU C (3) color singlet clusters of quarks and gluons contribute [TUIITTIIT^I 

I13j to the partition function of the system with their masses rrii and volumina 

Vi. 

The SBM with proper volume was solved analytically by the Laplace transform 
to the isobaric ensemble and, indeed, its solution showed various possibilities 
for the phase transition between the QGP and the hadron gas. Since then this 
technique has been successfully applied not only to solve far more sophisticated 
versions 111! 1121 1131 ITS] of the SBM, but also to find an analytical solution of a 
simplified version of the statistical multifragmentation model |16U17| for the nuclear 
liquid-gas phase transition. 

However, up to now the formulation of the SBM has some severe problems. The 
first one is the absence of a width for the heavy resonances. From the Particle Data 
Group we know that heavy resonances with masses to > 3.5 GeV may have 
width comparable with their masses. Taking the widths into account shall effectively 
reduce the statistical weight of the resonance. As we shall show below, this change 
may eventually remove the divergence of the SBM thermodynamics. The second 
problem arises while discussing the results of the hadron gas (HG) model [T91l2Tj] . 
The HG model describes remarkably well the light hadron multiplicities measured 
in nucleus- nucleus collisions at CERN SPS and BNL RHIC energies. This 
model is nothing else as the SBM of light hadrons which accounts for the proper 
volume of hadrons with masses below 2.5 GeV, but neglects the contribution of the 
exponentially growing mass spectrum. 

In other words, in order to calculate particle ratios within the HG model it is, on 
the one hand, necessary to consider all strong decays of resonances according their 
partial width collected in |18| . and, on the other hand, it is necessary to truncate 
the hadron spectrum for masses above 2.5 GeV. Thus, one immediately faces the 
following problem: "Why do the heavy resonances with masses above 2.5 GeV 
predicted by the SBM not appear in the particle spectra measured in heavy-ion 
collisions at SPS and RHIC energies?" Note that the absence of heavy resonance 
contributions in the particle ratios cannot be due to the statistical suppression of 
the Hagedorn mass spectrum because the latter should not be strong in the quark- 
hadron phase transition region, where those ratios are believed to be formed [T§ll2Tj] . 

In the present contribution we suggest that the introduction of a finite width 
of the resonances can solve the above problems of the SBM. In the next section 
we formulate a simple statistical model that incorporates besides of the Hagedorn 
mass spectrum also medium dependent resonance widths due to the hadronic Mott 
effect, and analyze its mathematical structure. In Section 3 we discuss a model fit 
to recent lattice data of QCD thermodynamics [TJ and some possible consequences 
for heavy- ion physics.. 
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2. Resonance Width Model: Mott Transition 



According to QCD, hadrons are not elementary, pointlike objects but rather color 
singlet bound states of quarks and gluons with a finite spatial extension of their 
wave function. While at low densities a hadron gas description can be sufficient, at 
high densities and temperatures, when hadronic wave functions overlap, nonvanish- 
ing quark exchange matrix elements between hadrons occur in order to fulfill the 
Pauli principle. This leads to a Mott- Anderson type derealization transition with 
frequent rearrangement processes of color strings (string- flip [23 ) so that hadronic 
resonances become off-shell with a finite, medium-dependent width. Such a Mott 
transition has been thoroughly discussed for light hadron systems in [SJJ and has 
been named soft deconfinement. The Mott transition for heavy mesons may serve as 
the physical mechanism behind the anomalous J/ip suppression phenomenon |23| . 

We introduce the width T of a resonance in the statistical model with the Hage- 
dorn mass spectrum through the spectral function 



A(s,m) = N s 



r m 



(s -m 2 ) 2 +r 2 



(1) 



a Breit-Wigner distribution of virtual masses with a maximum at y/s = m and the 
normalization factor 



N s = 



ds- 



Tm 



r 2 r 



i 



§ + arctan , 



(2) 



The energy density of this model with zero resonance proper volume for given 
temperature T and baryonic chemical potential \i can be cast in the form 

£=7T,p,. . . 

/OO pOO 
dm / ds p H (m) A(s,m) £a(T, ha] Vs), (3) 

where the energy density per degree of freedom with a mass m is 

d 3 k 



£A(T,HA;m) 



Vk 2 + m? 



(27r)3 exp (^&V) +( 5/ 



(4) 



with the degeneracy gA and the baryonic chemical potential (j>a of hadron A. For 
mesons, 5m = —1 , fiM = and for baryons 5b = 1 and /j,b = respectively. 
According to Eq. (jSJ) the energy density of hadrons consists of the contribution of 
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light hadrons for m,i < itia and the contribution of the Hagedorn mass spectrum 
Pnim) for m > itia- 

A new element of Eq. J3J in comparison to the SBM is the presence of the y/s- 
dependent spectral function. The analysis shows that, depending on the behavior 
of the resonance width T in the limit m — > oo, there are the following possibilities: 

• For vanishing resonance width, T = 0, Eq. Q evidently reproduces the usual 
SBM. 

• For final values of the resonance width, T = const, Eq. diverges for all 
temperatures T because, in contrast to the SBM, the statistical factor in Eq. 
© behaves as {exp [(mf, — /^)/T] + 5a} 1 so that it cannot suppress the 
exponential divergence of the Hagedorn mass spectrum pnirn). 

• For a resonance width growing with mass like the Hagedorn spectrum T ~ 
Cr exp ^ or faster, Eq. ® converges again. 

Indeed, in the latter case the Breit-Wigner spectral function behaves as 

Tm 



(s — m 2 ) 2 + r 2 m 2 



2 exp (-£-) (5) 



7r r \ t h j 



and cancels the exponential divergence of the Hagedorn mass spectrum. Hence, the 
energy density remains finite. Note that both the analytical properties of model 
@ and the right hand side of Eq. (J5J remain the same, if a Gaussian shape of the 
spectral function is chosen instead of the Breit-Wigner one. 

It can be shown that the behavior of the width at finite resonance masses is not 
essential for the convergence of the energy density J3J). In other words, for a con- 
vergent energy density above Th it is sufficient to have a very small probability 
density JSJ (or smaller) for a resonance of mass m to be found in the state with the 
virtual mass ^/s. Since there is no principal difference between the high and low 
mass resonances, we can use the same functional dependence of the width T for all 
masses. Thus, for the following model ansatz 



0, for T<T, 



H 



. ,N T ( 6 ) 



for T > T, 



H 



the energy density is finite for all temperatures and the divergence of the SBM 
is removed. At T — Th, depending on choice of parameters, it may have either a 
discontinuity or its partial T derivative may be discontinuous. As discussed above, 
for T < Th such a model corresponds to the usual SBM, but for high temperatures 
T > Th it remains finite for a wide choice of powers N m . 
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Figure 2: Fit of the lattice QCD data £Q with the Mott-Hagedorn resonance gas 
model ©■ For details see text. 

3. Applications for lattice QCD and heavy-ion collisions 

As one can see from Fig. ^the Hagedorn gas model correctly reproduces the 
lattice QCD results below the critical temperature T c and just in a vicinity above 
T c , but not for large temperatures. Fig. |3 shows a comparison of the same lattice 
QCD data pp with the Mott-Hagedorn gas 10 where the parameters of the spectral 
function are N T = 2.325, N m = 2.5 and T H = 165 MeV and m a = m b = 1 GeV. 
The successful description of the lattice energy density |T| indicates that above T c 
the strongly interacting matter may be well described in terms of strongly cor- 
related hadronic degrees of freedom. This result is based on the concept of soft 
deconfinement and provides an alternative to the conventional explanation of the 
deconfinement transition as the emergence of quasifree quarks and gluons. 

Another interesting feature of the model I© is that it allows to explain naturally 
the absence of heavy resonance contributions to the particle yields measured at 
highest SPS and all RHIC energies, where QGP conditions are expected ^OlED]- In 
order to find out whether a given resonance has a chance to survive till the freeze- 
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out it is necessary to compare its lifetime with the typical timescale in the system. 
There are two typical timescales usually discussed in nucleus-nucleus collisions, the 
equilibration time r eq and the formation time tj . The equilibration time tells when 
the matter created in collision process reaches a thermal equilibrium which allows 
to use the hydrodynamic and thermodynamic descriptions. For Au + Au collisions 
at RHIC energies it was estimated to be about r eq w 0.5 fm On the other hand 
in transport calculations the formation time is used: the time for constituent quarks 
to form a hadron. The formation time depends on the momentum and energy of 
the created hadron, but is of the same order r/ ~ 1 — 2 fm as the equilibration 
time. 

Since within our model the QGP is equivalent to a resonance gas with medium 
dependent widths, all hadronic resonances with life time r _1 (m) shorter than 
max{ry, r eq } will have no chance to be formed in the system. Therefore, the upper 
limit of the the integrals over the resonance mass m and over the virtual mass y/s 
in Eq. © should be reduced to a resonance mass defined by 

r(r7i) _1 = max{r/,T e(Z } . (7) 

This reduction may essentially weaken the energy density gap at the transition 
temperature or even make it vanish. Thus, the explicit time dependence should be 
introduced into the resonance width model © while applying it to nuclear colli- 
sions, and this finite time (size) effect, as we discussed, may change essentially the 
thermodynamics of the hadron resonances formed in the nucleus-nucleus collisions. 

4- Conclusions and Remarks 

The statistical bootstrap model allows to interpret the QGP as the hadron res- 
onance gas dominated by the state of infinite mass (and infinite volume). We argue 
that it is necessary to include the resonance width into the SBM in order to avoid 
the contradiction with the experimental data on hadron spectroscopy. We found 
that the simple model ©-© with a vanishing width below Hagedorn temperature 
Th and a Hagedorn spectrum-like width above Th may not only eliminate the di- 
vergence of the thermodynamic functions above T#, but it is able to successfully 
describe the lattice QCD data yQ for energy density with three fitting parameters 
only. Such a model also allows to naturally explain the absence of heavy resonance 
contributions in the fit of the experimentally measured particle ratios at SPS and 
RHIC energies. 

However, such a modification of the SBM requires an essential change in our 
view on QGP: it is conceivable that hadrons of very large masses which should be 
associated with a QGP cannot be formed in nucleus-nucleus collisions because of 
their very short lifetime. 

It is also necessary to remind that presented model should be applied to analyze 
the experimental data with care: it can be successfully applied to describe either 
the quantities associated with the chemical freeze-out, i.e. particle ratios or spectra 
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of f2 hyperons, <f>, J/tp and i/j' mesons that are freezing out at hadronization |26l 
E3E3ES|. But as discussed in Refs. [301 EH E21 the model presented here should 
not be used for the post freeze-out momentum spectra of other hadrons produced 
in the nucleus-nucleus collisions. Perhaps only such weakly interacting hadrons like 
f2, <f), J/ip and tp' will allow us to test the model presented. 

The question how to derive the non-zero width below Tjj has to be investigated. 
To solve this problem it will be necessary to include into present model a non- 
zero proper volume of hadrons. This, however, will require to consider a mixture 
of hadrons of different sizes as studied recently in [33[ EH or even a relativistic 
modification of the excluded volume of hadrons |341 135] . 
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